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ABSTRACT 

Synchrotron self-Compton (SSC) scattering is an important emission mechanism in 
many astronomical sources, such as gamma-ray bursts (GRBs) and active galactic 
nuclei ( AGNs) . We present a complete presentation of the analytical approximations 
for the Compton scattering of synchrotron emission with both weak and strong syn- 
chrotron self-absorption. All possible orders of the characteristic synchrotron spectral 
breaks (V a , u m , and v c ) are studied. In the weak self-absorption regime, i.e., v & < v c , 
the electron energy distribution is not modified by the self- absorption process. The 
shape of the SSC component broadly resembles that of the synchrotron component, 
but with the following features: (1) The SSC flux increases linearly with frequency 
up to the SSC break frequency corresponding to the self-absorption frequency i^ a ; (2) 
The presence of a logarithmic term in the high-frequency range of the SSC spectra 
makes it harder than the power-law approximation. In the strong absorption regime, 
i.e. u a > u c , heating of low energy electrons due to synchrotron absorption leads to 
pile-up of electrons to form a relativistic Maxwellian component besides the broken 
power-law component. According to such an electron spectrum, we calculate both the 
synchrotron and the SSC photon spectral components. We find that one rounded ther- 
mal bump appears in both spectral components. An analytical method is developed 
under the broken power law approximation, which is compared with the full numerical 
calculations. Due to the non-power-law behavior near the thermal bump, the analyt- 
ical approximation is not accurate around the transition region between the thermal 
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and non-thermal components. For the spectral regime v c < v & < f m , the spectrum is 
thermal-dominated (non-thermal-dominatcd) when v a > yju m v c {v a < ^Jv m v c ^) is sat- 
isfied. The analytical approximation for the SSC component is generally good enough 
at high energies above the non-thermal peak. For the thermal-dominated case, and 
the regime when v a is higher than both v m and v c , full numerical calculations are 
needed to reach an accurate result. We discuss possible applications and limitations 
of our approach to study astrophysical objects. In particular, we discuss an example 
of GRB reverse shock emission in a wind medium, which invokes v a > max(zy m , v c ). 

Key words: gamma ray bursts: general - radiation mechanisms: non-thermal 



1 INTRODUCTION 



Astrophysical sources powered by synchrotron emission should have a synchrotron self-Compton 
(SSC) scattering component. The same electrons that radiate synchrotron photons would scat- 
ter these synchrotron seed photons to high energies, forming a distinct spectral component. The 
SSC mechanism has been invoked to account for the observe d high energy emission in many 
astrophysical sources, such as gamma-ray bursts (GRBs) (e.g. 



1998 



Dermer et al. 



Wang et al. 



1998b 



2001 



2000 



Panaitescu fc Kumar 



Wu et a. 



Chiang fc Bottcher 



200. 



2002 



Zou et al.l 



2000 



Zhang et al. 



Meszaros et al. 



Sari fc Esin 



2001 



1994; 



Wei fc Lu 



Zhang fc Meszaros 



2009) and active galactic nuclei (AGNs) (e.g. 



2001 



Ghisellini et al 



201lh . 



SSC is a complex process. The flux at each observed frequency includes the contributions from 
electrons in a wide range of energy, which scatter seed photons in a wide range of frequency. 
Therefore, a precise description of the SSC spectrum invokes a complex convolution of the seed 
photon spectrum and electron energy distribution, which requires numerical calculations. However, 
for a synchrotron source with shock-accelerated electrons, the injected electron spectrum is usually 
assumed to be a simple power-law function, the corresponding electron energy distribution and 
seed synchrotron spectrum thus have simple patterns. Some analytical approximations for the SSC 
spectrum can be then made if the Compton scattering is in the Thomson regime. 

Besides the injected electron spectrum, two other factors are essential to define the shape of 
the final electron energy distribution in a synchrotron source: radiation cooling and self-absorption 
heating. There are three characteristic synchrotron frequencies in the spectrum: the minimum 
injection frequency (^ m ), the cooling frequency (z/ c ), and the self-absorption frequency (z/ a ). When 
u a < z/ c , the heating effect due to self- absorption is not important in modifying the electron 
energy spectrum. For a continuous injection of a power-law electron spectrum, the final electron 
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energy distribution is a broken power law. The seed synchrotron spectrum for SSC is characterized 
by a multi-segment broken power law, separated by v m , v c , and v a . Different ordering of the 
three characteristic frequencies leads to differe nt shapes of t he see d synchrotron spectrum. In the 
literature, usually v a < min(z/ m , v c ) is assumed. ISari fc Esinl (120011 ) have derived the approximated 
expressions of the SSC spectrum in the v & < u m < v c and z/ a < v c < u m regimes, respect iveljR 

When v & > v c happens, synchrotron self-absorption becomes an important heating source for 
the low-energy electrons. Consequently, the electrons are dominated by a thermal Maxwellian 
component until a "transition" Lorentz factor 7 t , above which the electrons are no longer a: 



by the self-a bsorption heating and keep the normal power law distribution (jGhisellini et al. 



1991 



are no longer a: 


'ected 


Ghisellini et al. 


19881 



1998al ). For these stron g absorption C8.SC S, ct rounded thermal peak would appear around v a 



in the synchrotron spectrum (IKobayashi et al. 
in the SSC spectrum. 



2004 ) , which would also result in some new features 



In this paper, we extend the analysis of ISari &: Esinl (120011 ) and present the full analytical 



approximated expressions of the SSC spectrum in all six possible cases of z/ a , z/ m , v c ordering. 
In Section 2, three weak synchrotron self-absorption cases (u a < u c ) are discussed. In Section 
3, we focus on the strong synchrotron self-absorption regime (z/ a > u c ), where synchrotron self- 
absorption significantly affects the electron energy distribution. We derive the expressions of both 
the synchrotron and the SSC spectral components, and compare the analytical approximations and 



the full numerical calculations. All the expres sions in th i s work are va ^ 
that the Klein-Nishina correction effect (e.g. 



Rees 



1967 



Nakar et al. 



id in the Thomson regime, so 



20091 ) is not important in the 



first order SSC component. We also li mit our treatment to the first-order SSC , and assume that 



the higher-order SSC components (e.g. 



Kobayashi et al. 



2007; 



Piran et al. 



20091 ) are suppressed by 



the Klein-Nishina effect. Such an assumption is usually valid for most problems. In order to make 
a simple analytical treatment, we have applied a simplified approximation for the synchrotron 
spectra, and adopted the simplification that the inverse Compto n scattering of m ono-energetic 



electrons off mono-energetic seed photons is also mono-energetic (ISari fc Esinl 120011 ). This would 



not significantly deteriorate precision of the analysis, while making it much simpler. 



1 Assuming weak self-absorption, Gou et al. (2007) derived analytical approximations of the SSC component for several other spectral 
regimes. 
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2 WEAK SYNCHROTRON SELF- ABSORPTION CASES 



In the single sc attering regime, the inverse Compton vo 



of electrons is ( IRybicki fc Lightman 



1979 



Sari fc Esin 



ume e missivity for a power-law distribution 



20011) 



J V 



3(77 



d 7 iV(7) / dxg(x)f Us (x), (1) 

7m JO 

where x = vj^ 2 v s (an angle-dependent parameter), f Us is the incident specific flux in the shock 
front, a? is Thomson scattering cross section, and g[x) = 1 + x + 2x In x — 2x 2 takes care of the 



angu 



ar dependence of the scattering cross section in the limit of 7 S> 1 (IBlumenthal fc Gould 



1970l ). One can approximate g(x) = 1 for < x < xq to simplify the integration, which would 
yield a correct behavior for 1 C 1 (Sari & Esin 2001). With such a simplification, the SSC spectrum 
is given by (Sari & Esin 2001), 



f 

j i' 



ic 



Ro 



T 



d 7 JV( 7 ) / dxf^x), 



(2) 



where (x) is the synchrotron flux, R is the co-moving size of the emission region, and the value 
of the parameter xq is set by ensuring energy conservation, i.e. x g(x)dx = J* x dx. 
When z/ a < u c , in the slow cooling regime (7 m < 7c ), the electron energy distribution is 

n (P ~ 1 )7ST 1 7" P ) 7m < 7 < 7c 



iV( 7 ) 



U (P- !)7 m X 7c7" 



p-i 



(3) 



7 > 7c- 



Here 7 m is the minimum Lorentz factor of the injected electrons, and p is electron spectral index. 
Cooling is efficient for electrons with Lorentz factor above the critical value 7c . Notice that EqE] 
is only valid for p > 1 . 

In the fast cooling regime ( 7c < 7m ), the electron energy distribution is 

2 



N(j) 



(4) 



«7c7 7c < 7 < 7m, 

^7m" 1 7c7" p_1 , 7 > 7m- 

In this regime, all the injected electrons are able to cool on the dynamical timescale. Therefore, 
there is a population of electrons with Lorentz factor below the injection minimum Lorentz factor 

7m- 

The seed synchrotron spectrum f Vs has spectral beaks at u a , u m and u c , where z/ a is the self- 
absorption frequency, below which the system becomes optically thick, and v m and v c are the 
characteristic sy nchrotron frequenc ies for the electrons with Lorentz factors 7 m and 7 C , respectively. 

As shown in lSari fc Esinl (l200lh . the critical frequencies in the SSC component are defined by 



different combination of 7 a ,7 m ,7c and z/ a , u m , u c . For convenience, we use a new notation in this 
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paper 

v \f = ^il v i x ^ i,j = a, c, m. (5) 

The physical meaning is the characteristic upscattered frequency for mone-energetic electrons with 
Lorentz factor 7$ scattering off mono-energetic photons with frequency Uj. 



2.1 Case I: z/ a < v m < v. 



This case has been studied by lSari fc Esinl (120011 ). The synchrotron spectrum read 



fu 



/max y — 
/max (j^ 
f max ( 7T" 



i-p 
2 



i-p 
2 



V > V r 



(6) 



where / max = f v {^m) is the peak flux density of the synchrotron component, which is taken as 
a constant. Subst ituting this seed photon spectrum into equation (j2J), the inner integral reads 
Isari fc EsiJbnOlh 



h — 9/max^o ( „^ 



V 



) /max ^0 



1 



47^y m a:o 



3 



(P+l) 



(P+2) 



47 2 f m a;o 



/max^O ( j, 



l-p 
2 



47 2 i/ c a;o 



< 47 2 z/ a x 
47 2 ^ a x < z/ < 47 2 ^ m x 
47 2 ^m^o < v < 47 2 z/ c x 
v > 47 2 z/ c x . 



(7) 



Similar to 



Sari fc Esinl ( 120011 ). only the leading order of v and zeroth order of ^a/^m and u m /u c 
are shown. However, we note that higher order small terms are needed to derive the following SSC 
spectrum (jSJ) through integrating the outer integral of equation (j2J). 

After integration, fl° is very complex. Keeping only the dominant terms, one gets the analytical 



2 Hereafter, the synchrotron spectra are denoted as fv(y) for simple presentation. Notice that when they are taken as seed spectrum, 
one should consider them as fv a {vs) and apply equation J2} to calculate the SSC spectra. 

© RAS, MNRAS 000, QUE] 



6 Gao, Lei, Wu, & Zhang 
approximation 

if ~ Ra T nf n 



x 



I max 2-0 

5 (P-1) f £a 

2 (p+1) V v « 

3 (p-1) 



(8) 



2(p-l/3) V^, 

(P-1) 
(P+1) 



"ma 

1 



IC . 



(p+1) 
(p- 



mm 



i-p 

2 



Notice that 



(p-i) f w 
(p+i) 



4(p+l/3) 
(p+l)(p-l/3) 

2(2p+3) 
(P+2) 



+ In 



Sari fc Esin 



2(2p+3) 
(P+2) 



(p+2) 
2 



+ In 



v < v, 



Z/ IC < V < z/ IC " 



Z/ IC < V < U lG ' 

^mc ^ " ^ ^cc i 



(P+2) 2 



(?+i)l n 

(P+2) 111 



IC 



( 1200 ll ) presented an opposite sign for the term 



(P+2) 5 



in the last segment, 



which might be a typo in that paper. 

The normalized synchrotron + SSC spectra for this and other two weak self-absorption cases 
are presented in Figure [TJ We note that these analytical expressions are not continuous around the 
breaks because of dropping the small order terms (see also Sari & Esin 2001), but the mis-match 
is small. When plotting the SSC curve in Figure [TJ we have used the analytical approximations, 
but added back some smaller order terms to remove the discontinuity. 



2.2 Case II: u m < u & < v c 



The synchrotron photons spectrum reads 



fn 
fn 
fn 
fn 



J.) + 4 
2 



2 



V 



1-P 

Vm 

1-P 
2 



v m < v < i/ a ; 
v a < v < z/ c ; 
v > u c ; 



(9) 



Evaluating the inner integral in equation we obtain 

47 2 z/ a x < f < 47 2 z/ c x 



1% — p^_i /max-^-O 



47 2 i/ m xo 
1 

(p+2)^ maxX ° Ut 



v 

47 2 !/ m xo ' 
l-p 

2 



(10) 



47 2 i/ m xo 



An interesting feature of this result is that I\ is linear with v all the way to v — 47 2 z/ a £ , 
indicating that a break corresponding to the break in the synchrotron spectrum at u m does not 
show up in the SSC spectrum for monoenergetic electron scattering. When v > 47 2 ^ a xo, the SSC 
spectrum follows the same frequency dependence as the corresponding seed synchrotron spectrum. 
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After second integration, we get the analytical approximation in this regime: 



ic 

v 



X 



Ro- T nf max x 



2(p+4)(p-l) f ^ 
3(p+l)2 



(pH 


-1) 


(p- 


-1) 


(pH 


-1) 


(p- 


-1) 


(P+l) 


(P- 


-1) 


(pH 


-2) 



mm 



mm 



P + l 

2 



1-P 
2 



1-P 
2 



1-P 
2 



1/ 

mm 



2(2p+5) 
(p+l)(p+4) 



+ In 



^ma 



2 + + 



2(2p+l) 



(P+l) 
_£ r- 

2 



+ ln p*- 



2(2p+5) , i , „ 
(p+2) + 111 1 -IC 



IC . 



U lC < V < Z/ IC " 

^ma ^ " ^ "mo 



I/ IC < 1/ < Z/ IC " 



^ca ^ " ^ ^cc ) 



V > V, 



IC 



Similar to the J result, there is no spectral break around v^m- Another comment is that the 
logarithmic terms make the SSC spectrum harder than the simple broken power-law approximation 
above the vF v peak frequency. At high frequencies, the simple broken power-law approximation 
may not be adequate to represent the true SSC spectrum. 



2.3 Case III: za < v r < v n 



This case was also studied by ISari fc Esinl (120011 ) . The seed synchrotron spectrum reads 



f»=< 



fmax\J,l) Ua 

/max (yj^ 
/max 

/max I 77"" ) ( ~ 



v <v a 



v c <v < z/ m ; 
, v > u m ; 



This gives 



J-l — 2 /max 3^0 ( „° 





2 /max^O ( 




3/max^O ( 




(p+2) / max 



47 2 J/ C 2!0 



47 2 i/ c x 



1 



v < 47 2 i/ a x 
47 2 ^ a a;o < ^ < 47 2 z/ c x 
47 2 i/ c x < ^ < 47 2 z/ m x 
z/ > 47 2 zv m x 
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and the final SSC spectrum 



f 



ic 



x 



Ra T nf n 

5 



_9_ 

10 

1 
3 

1 
3 



v 



7JC" 



AC 



1 



(p+2) 1% 



2(p+5) 
(p+2) (p-1) 
p 



,.ic 



2(P-1) 
3(p+2) 



hi 



JC 
m m 



2 (p+5) _ 2 (p-1) 

3 (p-1) 3 (p+2) 



+ In 



f < 


^ca ) 




-ca c 


< V < 


IC. 

^CC ' 




< V < Z/ IC " 


^IC 

cm 


< Z/ < 


^IC . 

^mmi 


V > 


z/ IC . 

mm 





(14) 



We note that ISari fc Esinl (120011 ) has an opposite sign in the term In y^hj i n the third segment, 
which might be another typo in that paper. 

We define ratio between the SSC luminosity and the synchrotron luminosity as the X parameter 
similar to Sari & Esin (2001), i.e., 



X 



'IC 



'syn 



P h 



(15) 



where {7 p h and Ub are the synchrotron photon energy density and magnetic field energy density, 
respectively. 

For i/ a < v m < u c (case I) and u m < i/ a < v c (case II), the vf v peaks of the synchrotron and the 
SSC components are at v c and respectively (see Figure 1). One can estimate 



X 



'IC 



'syn 



cc J v \ cc / 

ic \ — 

\ "mm J 



"cjmax I ,. 



1-p 

2 



1-p 



(16) 



which is consistent with Sari & Esin (2001). Note that when calculating X, we did not include the 
coefficients in the analytical approximations of the SSC component, which is of order unity. 

For z/ a < u c < u m (case III), the vf v peaks of the synchrotron and SSC components are at z/ m , 
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2 





2 




log 10 v[Hz] 

Figure 1. Total synchrotron + SSC spectra for weak synchrotron reabsorption cases (y a < u c ). The top panel is for i/ a < v m < 
case; the middle panel is for v m < i/ a < "c case; and the bottom panel is for f a < v c < case. The thin solid line is synchrotron 
component. The thick solid line in the SSC component is drawn using the analytical approximations, while the dashed lines are the 
broken power-law approximation for comparison. In all the cases, the vF v peaks for both the synchrotron and the SSC components are 
normalized to unity. 
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and v ^m i respectively. One therefore has 



X " IC 



,IC flC^IC ) 
mmJ v \ mm) 



^m.J 'max 

,2 



Cm 



~ 4x a T n J R7 c 7 m , (17) 
which is also consistent with Sari & Esin (2001). 



3 STRONG SYNCHROTRON SELF-ABSORPTION CASES 

When z/ a > z/ c , synchrotron/SSC cooling and self- absorption heating would reach a balance around 
a certain electron energy. The electron energy distribution would include two components: a ther- 
mal component and a non-thermal power-law component. We introduce a relativistic Maxwellian 
distribution with a peak at 7 a to denote the thermal component, which is superposed on the broken 
power law component. The power- law indices of the non-thermal part depend on the order of z/ a , 
u m , and v c . 

For the case of v c < v & < u m , the electron distribution becomes 



n^re •» 7c < 7 < 7t , 



^7c7~ 2 > 7t < 7 < 7m, ( 18 ) 

^7m" 1 7c7" p_1 , 7 > 7m- 

where the "transition" Lorentz factor j t can be solved numerically based on conservation of the 
electron number. 

For the case of u m < u c < z/ a , one has 

{ 4 7 2 -22 
7a, 7m<7<7t, , . 

W (P- !)7m 7c7 p , 7 > 7t- 
For the case of v c < u m < z/ a , one has 

n 7 a 3 6 7a ' 7c_7_7t, ^ 
^7m" 1 7c7" p " 1 , 7 > It- 

Following these new shapes of the electron distribution, the synchrotron photon spectra can 
be calculated, which also contain a thermal component and a (broken) power-law component. 
Still applying equation ©, one can both numerically and analytically calculate the SSC spectral 
component for another three cases in this regime. 
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3.1 Case IV: v c < v & < v m 

In this case, the synchrotron photon spectrum reads 



f f ( 

J max I 



''a 



2 



3 1 

8 , _1 (21) 



mi 



where a rounded thermal peak appears around z/ a (IKobayashi et al.l I2004T 

For the Maxwellian prescription of the electron and photon distributions, the SSC component 
for this regime can be only obtained numerically using equation ([2]). In order to derive an analytical 
approximation, one needs to make further simplifications. We attempt to simplify the electron 
distribution and synchrotron photon spectrum with the following broken power-law form: 



Niri) 



n-3r, 7 < 7a, 

'a 



^7c7 2 , 7a < 7 < 7m- (22) 



"7m X 7c7 p 1 , 7 > 7m- 



and 

2 



/max (^) , V < I/ a ; 

/maxf(^)^, ^a<^<^m; (23) 

f » (WT* Z y> Z y ■ 

jmax 4 , ^"mi 

where £H is the discontinuity ratio in the electron distribution at 7 a , 



fJK = — . (24) 

7a 

With such an approximation, one can derive 

' h ^ U ax x + 1) (jpfc) , * < 47 2 ^xo 

ia - l/max^o^H (4^^) 2 , 4 7 2 i/ a x < 1/ < 4 7 2 i/ m x (25) 

1 _ p 

J 3 - 2(^2) /maxXQ^ {jA * (4^— ) ' , U > 4rfv m X 



I = < 



3 The numerical value 3 in the exponent of the exponential term is designated to assure that the peak flux of the thermal component 
to be at v a . 
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and 



f, 



ic 



X 



Ro-Tnf mSkX x 
' !)(!* + 4) 



21 

■i 



7JU 



v 
7JC 



~ 5 



77TCT 



TjcT 



-2- - I + I In 

p-l 3 T 2 111 



1 m 2 I i^s, 
2(p+2)-' 1 l;y m 



mm 



,ic 

m m 



4 
p+3 



p-l 



Til 

7c 



3(P+1) 
(p-l)(p+2) 



(26) 



v < v. 



IC. 



I/ IC < 1/ < Z/ IC " 

^aa ^ " ^ ^am) 



^am ^ " ^ ^mm> 



IC 



In this case, there are two peaks in the ^i 7 !, spectrum for both the synchrotron and the SSC 
components. For the synchrotron component, the thermal peak is at (25/9)z/ a ~ 2.8u a , and the 
non-thermal peak is at u m . For the SSC component, the thermal peak at , and the non-thermal 
peak at v^m- The relative importance of the two peaks depend on the relative location of is a with 
respect to u c and u m . More specifically, the spectrum is non-thermal-dominated when z/ a < y 7 !^!^, 
and is thermal-dominated when z/ a > ^v ni v c - 

In Figure 2, we compare the full numerical result (solid) of the SSC component with the above 
simplified analytical approximation (dashed) as well as a simplest power law prescription (dot- 
dashed) of the SSC component. The non-thermal-dominated and the thermal-dominated cases 
are presented in Fig. 2a and Fig. 2b, respectively. Below v^L, The analytical approximation is 
obviously lower than the numerical one. This is due to the poor broken power-law approximation 
for the thermal bump. Above the analytical approximation fits well with the numerical result. 
Similar to the weak self-absorption regime (cases I-III) , the logarithmic terms make the analytical 
spectrum harder than the simple broken power-law approximation above the non-thermal vF v 
peak frequency. At high frequencies, the simple broken power-law approximation is not adequate to 
represent the true SSC spectrum. As shown in Fig. 2b, in the thermal-dominated case, the analytical 
approximation becomes poor. In this regime, one should use the full numerical calculation to obtain 
an accurate result. 
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3.2 Case V and VI: v & > max(z/ m , z/ c ) 

For these two cases, the treatments and results are rather similar to each other, we take y m < v z < 
z/ a as an example. In this case, the synchrotron spectrum reads 



/max (y u& 



/max - e J . 



1 3 -3f^)^ 



/n 



z/ a < V < U t ] 



v > u t ; 



(27) 



Applying equation one can numerically calculate the SSC component. 

Nonetheless, similar to case IV, we attempt to simplify the electron distribution and syn- 
chrotron photon spectrum as 



la. 



n ip ~ !)7 m X 7c7" 



/max y Va 



N(rf) 

and 

fu = 

where 



7a \7a 



p-1 



7 < 7a, 
7 > 7a- 



(28) 



JmaxT I 7T ) ' ^ ^ ^a 



(29) 



p-i 



(30) 



The inner integral 7 can be then approximated as 

f h ^ /maxXQ (j^jj + A > " < 

1 — \ v _£ 



^ 2 — 2fm+2'i /max^O^ 



2(p+2) ^ max-u- * ^4r i i/»x Q 

Integrating over the outer integral, one gets 



(31) 



IC 



Ro- T nf m3uX x 



(32) 



x 



2(p+2) 



1 

p+2 



>(^ + 4 



7JU , 



p+3 ~ ^ l V 2 (P+2) 



V > V. 



IC. 

aa > 

IC. 

aa ) 



Comparing the numerical result with the analytical approximation (Fig. 2c), we find that the 
analytical approximation cannot catch the key feature of the spectrum. This is because in the 
regime, the spectrum is thermally dominated. We suggest that for this regime, one should use full 
numerical calculations to obtain an accurate result. 

The case of v c < v m < u a is almost identical to the above v m < v c < u a . The only difference is 
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that the expression of %\ is modified to 
p-i 

- 

7a \7a/ 

This is understandable, since in the fast cooling, the electron energy spectral index is p = 2, so 
that the factor (p — 1) can be reduced to 1. The numerical and analytical results in this regime is 
identical to Fig.3c. Again one needs full numerical calculations to obtain accurate results. 
Finally, we investigate the X parameter in the strong synchrotron self-absorption regime. 
For v c < v a < v m (case IV), if the spectrum is non-thermal-dominated, the synchrotron and 
SSC emission components peak at u m and respectively. One thus has 
t ,.ic fiCY,.ic \ 



X 



v^ m Ra T nf max x W- 



i 

,IC \~2 

mm 



~ 4xla T nR-f c -f m . (34) 

If the spectrum is thermal-dominated, the synchrotron and SSC emission components peak at z/ a 
and respectively. One has 

t ,.IC flC(,.lC\ 



X 



~ 4x 2 a T nRjl (35) 

In general, the X parameter for v c < v a < u m (case IV) is Ax^arnR ■ max(7^, 7 c 7 m )- 

For u m < v c < u a (case V) and v c < u m < v a (case VI), the synchrotron and SSC emission 
components peak at v a and u^, respectively. In this case, one has 

t ,ICflC(„IC\ 



X 



^Ra T nf m ^x 



^a/max 

~ AxlamR^. (36) 

This is the same as the thermal-dominated case for u c < z/ a < z/ m (case IV), we thus conclude that 
the expression of X is equation (136)) as long as the spectrum is thermal-dominated. 

© RAS, MNRAS 000,[TJ{T7] 



Compton Scattering of Self- Absorbed Synchrotron Emission 15 




Figure 2. Same as Figure [T] but for strong synchrotron reabsorption cases . The top panel is for u c < < u m and < y/u m u c case; 
the middle panel is for u c < < z/ m and i/ a > \Ji>mPc case; and the bottom panel is for v & > max(i/ m , ^ c ) case. All the solid lines are 
numerical results, while the dashed lines are analytical approximations and the dot-dashed lines are broken power-law approximations. 
© RAS, MNRAS OOO.HHIZ] 
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4 CONCLUSION AND DISCUSSION 



We have extended the analysis of ISari fe Esinl (120011 ) and derived the analytical approximations 
of the SSC spectra of all possible orders of the three synchrotron characteristic frequencies v a , is m , 
and u c . Based on the relative order between u a and u c , we divide the six possible orders into two 
regimes. 

In the weak self- absorption regime z/ a < u c , self-absorptio n does not af f ect th e electron energy 
distribution. Two cases in this regime have been studied by 



Sari fc Esin 



( 20011 ). Our results are 



consistent with theirs (except the two typos in their paper). For the other regime u m < z/ a < u c , we 
find that the SSC spectrum is linear to v all the way to u^ a , and there is no break corresponding 
to z/ c m . 

In the strong self-absorption z/ a > v c regime, synchrotron self-absorption heating balances syn- 
chrotron/SSC cooling, leading to pile-up of electrons at a certain energy, so that the electron energy 
distribution is significantly altered, with an additional relativistic Maxwellian component besides 
the non-thermal power law component. Both the synchrotron and the SSC spectral components 
become two-hump shaped. To get an analytical approximation of the SSC spectrum, we further 
simplified the electron energy spectrum and the seed synchrotron spectrum into a broken power- 
law form, with a sharp cutoff above the piling up energy. For the v c < z/ a < v m case, the analytical 
approximation fits well with the numerical result in high energy range above the non-thermal 
vF v peak, but is a poor representation of the numerical results near the transition region from 
the thermal to the non-thermal component. We suggest that for the thermal-dominated cases, i.e. 
^ a > v^m^c in the z/ c < v, A < v m regime or the v a > max(z/ m , u c ) regime, full numerical calculations 
are needed to get accurate results. 

In general, the SSC component roughly tracks the shape of the seed synchrotron component, but 
is smoother and harder at high energies. For all the cases, we compare our analytical approximation 
results of SSC component with the simplest broken power-law prescription. We find that in general 
the presence of the logarithmic terms in the high energy range makes the SSC spectrum harder than 
the broken power-law approximation. One should consider these terms when studying high energy 
emission. The only exceptions are the z/ a > max(V m , v c ) regimes. However, in these regimes the 
analytical approximations are no longer good, and one should appeal to full numerical calculations. 

Our newly derived spectral regimes may find applications in astrophysical objects with high 
"compactness" (i.e. high luminosity, and small size). In these cases, v a can be higher than v c or u m , 
or even both. For example, in the early afterglow phase of GRBs, slow cooling may be relevant, and 
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the radio afterglow is self-absorbed with z/ a above u m (e.g. 



Chandra fc Frail 



20121 ). In the prompt 



emis sion phase when fast cooling is more relevant, the self-absorption frequency can be above u c 



(e.g. 



Shen fc ZhaneJ 120091 ) . 



An example of the extreme case z/ a > max(z/ m , u c ) can be identified for a GRB problem. For a 
dense circumburst medium with a wind-like (n cx r~ 2 ) structure, in th e reverse shock region, the 
condition u a > max(i/ m , u c ) can be satisfied (e.g. 



Kobayashi et al. 



20041 ). For a GRB with isotropic 



energy E = 10 52 E 52 , initial Lorentz factor T = 100r 0j2 , initial shell width A = 10 12 A 12 running 
into stellar wind with density p = (5 x 10 11 g cm _1 )A*r -2 , one can derive following parameters 
at the shock crossing radius r x : The blastwave Lorentz factor T x = 25.8A* 1 ^ 4: A^^eU 4: , u, 



3.1 x 10 14 Hz [g(p)/g(2.3)]A,A- 2 1/2 E l 



o2 



1/2 p 2 
e e,-l e B,-2 i 0,2' V c 



1.2 x 10 12 Hz A- 2 A!/ 2 ^ /2 e- 3/2 



12 ^52 5 "m 
12 -^52 c R-2' 



4.6 x 10 14 Hz v4: /J A 



3/5 A -11/10^,1/10^3/10 r -2/5 



12 



52 



R-2 1 0,2 



Here e e = 0.1e e _i and e# = 0.01e#_2 are microphysics 



shock parameters for the internal energy fraction that goes to electrons and magnetic fields, p is 
the electron spectral index, and g(p) = (p — 2)/(p — 1). We can see that for typical parameters, 
v a > max(j/ c , u m ) is satisfied. In this regime, one should check whether the "Razin" plasma effect 
is important. At shock crossing time, the comoving number density of the shocked ejecta region 
is n' = 2.3 x 10 8 cm" " 3 Al /4 A^ /A E; 2 1/4 Tol Noticing that the comoving plasma angular frequency 
is u' p = 5.63 x 10 4 s -1 ?^' 1 / 2 , one can write the plasma frequency in the observer frame as u p = 

102Ay 20 A- 1/20 ^2 /20 e ^/2 roi /5 , one gets 



1.4 x 10 11 Hz aTa-^E^T-^ 2 . Multiplying by 7a 



J h2 



-B.-2 1 0,2 ) 



1.4 x 10 13 Hz AV'^A 



17/40 a —47/40 j-i7/40 -l/10-p-7/10 



12 



-^52 ^3,-2 1 0,2 



, which is much smaller than v a . This suggests 



that the Razin effect is not important (jRybicki &: Lightmanlll979l ). and the dominant mechanism 
to suppress synchrotron emission at low energies is synchrotron self-absorption. 

Finally, the formalism of this paper only applies to the one-zone SSC models in the Thomson 
regime. The theory needs to be extended in order to study more complicated systems. 
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